arXiv: 1503.06907vl [hep-ph] 24 Mar 2015 


Eur. Phys. J. C manuscript No. 

(will be inserted by the editor) 


Study of Ah Al^l and A^ plu decays in the Bethe-Salpeter 
equation approach 


Y. LirpE, L.-L. LirPU, X.-H. GucpE 

^College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, People’s Republic of China 


Received: date / Accepted: date 


Abstract In our previous work, based on the SU{6) 
spin-flavor wave function, we regard A and p as com¬ 
posed of different quark-diquark configurations and es¬ 
tablished the Bethe-Salpeter (BS) equations of config¬ 
urations for quark and scalar diquark. In our present 
work, we apply this model to calculate the form factors 
of the semileptonic transitions At —>■ Al~^l~ (/ = /i, e, r) 
and At —)► pW within the Standard Model (SM). The 
decay At —>■ Ap^p~ is especially interesting since it has 
been measured in CDF and LHCb Collaborations and 
this rare decay is very sensitive to new physics effects. 
The decay At ^ pW is a, promising mode for the mea¬ 
surement of the Cabibbo-Kobayashi-Maskawa matrix 
element \Vub\ at the Large Hadron Collider. In our cal¬ 
culations, depending on the ranges of the parameters in 
the model including the diquark mass and the interac¬ 
tion strength between the quark and the diquark in the 
kernel of the BS equation, we find that the branching 
ratio of At —> ApApT in our model is consistent with 
the experimental data and the current experimental re¬ 
sults from LHCb agree with the differential branching 
ratio of At —>■ Ap'^p~ from our calculation except at the 
lager momentum transfer region. This indicates that 
there is still room for possible new physics effects. We 
also give comparisions of the total branching ratios of 
At —>■ AAl~ and At —?► plv with those given by other 
phenomenological methods. 


I. Introduction 

In recent years, a lot of experimental progresses have 
been made in spectroscopy and decays of heavy baryons 
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containing a heavy bottom or charm quark 
Processes associated with the flavor-changing neutral 
current 6 —> s transition have regained much atten¬ 
tion since the CLEO measurement of the radiative b —> 
S 7 decay [5]. Although the experimental measurement 
of mesonic b ^ s transitions appeared about twenty 
years ago [7] , the first observation of the baryonic decay 
At Ap'^p~ was reported in 2011 by the CDF Collab¬ 
oration [5]. A first measurement of the differential and 
total branching fractions for this rare decay by LHCb 
was reported in 2013 [H]. The decay At —?► AAl~{l = 
p,e,T) proceeds through electroweak loop diagrams in 
the Standard Model (SM). Since non-Standard Model 
particles such as supersymmetric particles m and light 
dark matter particals m may also participate in these 
loops, measurement of this decay can be used to search 
for new physics m- Furthermore, this channel can be 
used as a tool in the exact determination of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements, Vtt and 
Vts, and in the study of CP and T violations. On the 
other hand, the study of the exclusive decay At —> pW 
is a promising mode for the measurement of the poorly 
known magnitude of the CKM matrix element \Vut\ at 
the Large Hadron Collider. So far, all the measurements 
of |I4h| have been from B meson decays and were per¬ 
formed at B factories m- 

Theoretically, there are some works devoted to the 
analysis of At —?► AAl~ decays in the SM and in var¬ 
ious scenarios of physics beyond the SM [BIISIIISIIIZI 
[nnmi. In order to use the decay to search for new 
physics, one should determine the hadronic matrix ele¬ 
ment At ^ A, which is expressed in twelve form factors. 
In the heavy quark limit, with the application of the 
heavy quark effective theory (HEQT) for the b quark. 
At —>■ A(j)) transition can be described by two inde¬ 
pendent form factors |21j . The decrease in the num- 
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ber of form factors greatly simplifies calculations. How¬ 
ever, these two form factors contain all soft QCD effects 
which are difficult to calculate from the first principles. 
Therefore, one needs to resort to some phenomenologi¬ 
cal models. The ^ A form factors were calculated in 
quark models P^[^[Mll251l26l[?7] and the perturbative 
QCD approach [55]. Moreover, the simple pole model 
was adopted to compute the form factors by Manuel 
and Rochsiegel [53]. The authors of Refs. [15-16,29- 
37] employed the widely applied approach of QCD sum 
rules to calculate these two form factors. The informa¬ 
tion on Ac ^ A form factors is available from experi¬ 
mental measurement of the semileptonic Ac —^ Ae'^Ve 
decay [381139] . and this information was used to con¬ 
strain the Ab ^ A form factors Recently, the 

form factors of the At ^ A transition were determi¬ 
nated in the first lattice QCD simulation [JO]. In our 
previous work, in the “quark-diquark” model in which 
a baryon is regarded as a bound state of a quark and a 
diquark, we established the Bethe-Salpeter (BS) equa¬ 
tions for the quark and scalar diquark configurations 
of A and p based on SU (6) spin-flavor wave functions 
[msa. Then we solved them in the covariant instan¬ 
taneous approximation with the kernel containing both 
the scalar confinement and one-gluon-exchange terms 
[iTlB51H51l44ll451li5] . In the present work, we will apply 
the BS wave functions of A and p and those of Ai, which 
were obtained previously to calculate branching ratios 
for Ab —>■ Al~^l~ and Ab —>■ pli> and compare our results 
with experiment, lattice data, and results from other 
phenomenological methods. 

The layout of the paper is as follows, in Section II, 
we will review the basic formalism for the BS equations 
for A and p. With the aid of HQET, the form factors 
which are involved in Ab —>■ Al^l~ and Ab —> plv de¬ 
cays will be calculated using the BS wave functions of 
H, p and those of Ab which were obtained previously. In 
Section III we will show calculations and the numerical 
results of the branching ratios of Ab —>■ Al~^l~. In Sec¬ 
tion IV, we will give the branching ratios of Ab —>■ pW, 
and will compare with the results from QCD light-cone 
sum rule. Section V will be deserved for a summary and 
some discussions. 


II. Form Factors for Ab —>■ Al~^l and Ab —>■ pli^ 

A(jp) contains three light quarks u, d, and s (u), in which 
all the three light quarks play important roles in the 
dynamics inside the baryon. In general, the parity of a 
baryon at the ground state is positive. Since the parity 
of the quark in the baryon is supposed to be positive, 
the parity of the diquark involved in the ground state 
baryon should also be positive. Due to Pauli principle. 


two quarks with the same flavor constitutes an axial- 
vector diquark and two quarks with different flavors can 
constitute either a scalar diquark or an axial-vector di¬ 
quark. Based on the SU (6) wave function of the proton, 
the proton state can be expanded in the terms of quark- 
diquark configurations as follows m 




3^2 


[3u^{ud)ofi + u^{ud)ifi — 


V2u^{ud)i^i — V2d^{ uu)i^0 + 2d^{uu)ip], 


pi = 


3y/2 


[3u'*'(u(i)o,o ~ u^{ud)i^-i 


V2u^{ud)i^-i -|- '/2d^{ uu)l,0 — 2d^{uu)i^-i]. 


( 1 ) 

( 2 ) 


In the same way, we can obtain the following forms for 
A [H]: 

Hi = —p[2sI(itd)Q,o + •\/2(iI(Ms)ip — dl(us)ip + 

2-^3 

d^(us)o,o - '/2ul(ds)ip -I- ul(ds)ip - 
uI(ds)o,o] (3) 

HI = ^^[251(14(1)0,0 — '/2d^(us)i^-i + dI(iis)i,o 

-|-dI(iis)o,o + '/2u'^(ds)i^-i - Ml(ds)i,o - 
M'‘'(ds)o,o]- (4) 

In Eqs. (1-4) the first and the second subscripts corre¬ 
spond to the total spin and the third component of the 
spin of the diqaurk, respectively. The arrow f (1) indi¬ 
cates the spin direction of the corresponding baryon is 
up (down). 

Ab is regarded as a bound state of a 6 quark and 
a scalar diquark {ud)op. In order to calculate the form 
factors in the transition Ab ^ A (p), where the b quark 
decays into s (u) quark and the scalar diquark behaves 
as a spectator, one should calculate the BS wave func¬ 
tions of the configurations s(ud)o,o and ii(ud)o,o in A 
and p, respectively. 

We define the BS wave function of the system q{ud)ofi 
{q = s or u) as the following: 

x{xi,X2,P) = {Q\T'ii}{xi)ip{x2)\P), (5) 


where ipixi) and ifi{x 2 ) are the held operators of the 
light quark at position xi and the light scalar diquark 
at position X 2 ^ respectively, P = Mv is the momentum 
of A or p, and M (u) is its mass (velocity). Let 
and rriD represent the masses of the light quark and 
the light diquark in the baryon A or p, Ai = ^ ’ 

A 2 = , and p represent the relative momentum 

of the two constituents. X = AiXi -I- A 2 a :2 is the coor¬ 
dinate of the center of mass and x = Xi — X 2 - The BS 
wave function in momentum space, Xp(p), is related to 
x{xi,X 2 , P) through the following equation. 
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X{XUX2,P) = j xp{p)- ( 6 ) 

The BS equation for the q{ud)o,o system in momen¬ 
tum space can be written as follows: 


Xp{p) = Sf{pi) 


dU 


K{P,p,q) xp{q) Sd{p2), 


(7) 


where pi = XiP -I- p and p 2 = X 2 P — P are the mo¬ 
menta of the light quark q and the light scalar diquark, 
respectively, K{P, p, q) is the kernel which is defined 
as the sum of two particle irreducible diagrams, Spipi) 
and Sd{p 2 ) are propagators of the light quark with mo¬ 
mentum Pi and the light diquark with momentum p 2 . 
Motivated by the potential model, the kernel is given 
by miH] 


-iK{P,p,q) = I IVi{p,q) + Xtj. r^V 2 {p,q), (8) 


where P^ = {p 2 -b q 2 )^F{Q^) is the effective vertex of a 
gluon with two scalar diquarks, F{Q‘^) is introduced to 
describe the structure of the diquark [4^148] . F{Q‘^) = 
^ where Qq is a parameter which freezes F{Q'^) 
when is very small. In the high energy region the 
form factor is proportional to ^ which is consistent 
with perturbative QCD calculations. By analyzing the 
electromagnetic form factor for the proton, it was found 
that Qq = 3.2GeV^ can lead to consistent results with 
the experimental data. Vi and V 2 are the scalar con¬ 
finement and one-gluon-exchange terms which have the 
following forms in the covariant instantaneous approx¬ 
imation respectively 


Vi{pt 


-qt) 


Sttk 

[{pt - qtY + tPf 


-{2Trfd^{pt - qt) 


d^k Sttk 


(9) 


V 2 {pt - qt) = 


IGtt 

3 {pt - qt)'^ + ’ 


( 10 ) 


where pt and qt are the transverse projection of the 
relative momentum along the momentum P, which are 
defined as p^ = p^ — v ■ pv^ and q^ = q^ — v ■ qv^. The 
second term of Vi is introduced to remove the infrared 
singularity at the point pt = qt, and the small parame¬ 
ter p is introduced to avoid the divergence in numerical 
calculations. After considering the constraints on xpip) 
imposed by parity and Lorentz transformations, xp{p) 
can be expressed in terms of two Lorentz-scalar func¬ 
tions, /i and / 2 , 


( 11 ) 


where u(u, s) is the Dirac spinor of A or p. Defining 
/i( 2 )(= / ^/i( 2 )), we find these two BS scalar wave 
functions satisfy the coupled integral equations as fol¬ 
lows: 

I ^ _ I _ 

^ J (27r)3 4 w,w_d(M - Wq - wd) 

[{mg + ujq)(yi + 2uJDV2F{Q^))fi{qt) 

-{Pt ■ qt + Pt)V 2 F{Q'^)fi{qt)] 

_ I" dApt _1_ 

J (27r)3 4 w,w_d(M - Wq - wd) 

[iVi-2cuDV2FiQ^))pfqtf2iqt) 

-{mg F UJg){pt ■ qt + qt)V 2 F{Q^)f 2 {qt )], (12) 


MPt) = - 


d^Pt 


(27r)3 AujqUj£){M — LOq — iOp) 

[{Vi+2u:DV2F{Q^))fi{qt) 

-{rUg - ujg) ^^* V 2 T^(Q^)/i(gt)] 

Pt 

_ f dApt _1_ 

J (27r)3 4 w,w_d(M - Wq - wd) 

[{mg - u;g)(Vi - 2uJDV2F{Q^))^^f2{qt) 

Pt 


-{pt • qt + qt)y 2 F{Q'^)f 2 {qt)], 


(13) 


where ujd = \Jrn^ + Pt and Wg = y + pj. The 

BS wave functions /i( 2 ) can be solved numerically by 
discretizing the integration region (0, 00 ) into n pieces 
(n is chosen to be sufficiently large). The normalization 
condition for the BS wave function is given in the fol¬ 
lowing after imposing the covariant instantaneous ap¬ 
proximation on the kernel [4311^1^ : 




d^qd^p 

(27r)8 


xp{p, 


s) 


= Sss’xp{q, s') 


(14) 


where fi( 2 ) and ji( 2 ) represent the color indices of 
the light diquark and the light quark, respectively, s*- ^ 
is the spin index for the light baryon. 


iij2 ji h h 31 ^ 


/p* 2 ^ 2 ji jg inverse of the four point propagator 
defined as follows: 


q) = - q)sy^'>{pi)sy\p 2 ). 

(15) 

In our calculations, we choose the diquark mass mp 
to range from 700 MeV to 800 MeV [57]. With this 


Xp{p) = ifi +l/tf 2 )u{v,s), 
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choice for niD, the binding energy E is negative and 
varies from -90 MeV to -190MeV. k is chosen to range 
from 0.02 GeV^ to 0.08 GeV^ [27ll4^ . Then, for each 
mu, we get a value of aseff corresponding to a value of 
K. Solving the discretized Eqs. (12,13), which become 
an eigenvalue equation, we obtain the numerical results 
for /i(pt) and / 2 (pt), which depend on two parameters, 
mo and k. 

Using Lorentz symmetry and discrete C, P, T sym¬ 
metries, one can show that the following matrix ele¬ 
ments of the A}) ^ A transition can be parametrized by 
twelve independent form factors [35], 


{A{P',s')\-s^^h\A^{P,s)) = 

ua[P', + ig 2 (Jij.uP'' + 9 ^Pii)ua^{P, s ), 

{A{P',s')\s-ffj,-f5b\Ab{P,s)) = 

ua{P', + it 2 (Tt,^p'' + hp^)^bUAt{P, s), 

{A{P',s')\Maf^’'g''b\Ab{P,s)) = 

ua{p', s')(5'i’7^ + igl + gIq^)uAtiP. s), 

{A{p', s')\Maf^''j 5 g’'b\Ab{P, s)) = 

UA{P',s'){t'^-f^^ + + t'^qf^)^5UA,{P,s), (16) 

where q = P — P' is the momentum transfer, and gi, ti, 
gf ,tf (i = 1, 2 and 3) are various form factors which are 
Lorentz scalar functions of q^. The most general form 
for the matrix elements in Eq. (16) consistent with the 
spin symmetry on the b quark in the limit mb —^ oo is 


{A{P\s')\sr^b\ Ab{v,s)) = UAiP\s'){Fi + F 2 '^) 

r>^UA,{v,s), (17) 

where /), represent 7 ^, 7 ^ 75 , and cr^^^bq''. Fi 

(i = 1,2) can be expressed as functions solely of v ■ 
P', which is the energy of the A baryon in the Ab rest 
frame. Gomparing Eq. (16) with Eq. (17), we obtain 
the following relations: 


The transition matrix for Ab ^ A can be expressed in 
terms of the BS wave function of Ab and the s(ud)o,o 
component of Tl, Xp'ip')^ 


{AiP',s')\sr,\Ab{P,s)) =/^x 

xi'{p')rpXp^ip)Su^{p2)- (19) 

Erom Eqs. (3) and (4) we can see that the Glebsh- 
Gordan coefficient of the s(ud)o,o configuration is 1 /\/3. 
Substituting the BS wave functions of Ab and the s(ud)o,o 
system into Eq. (19) and integrating out the longitudi¬ 
nal momentum pi , we obtain the following forms for Fi 
and F 2 : 


Fi 


1 

7! 



2 w((Myi — (jjujp) — Lo'g — bfbP 


1 cos 9pt) 


{[(7(p) - kt) + 2uJD%(.p't - kt)F{Q^)){ms -h w') 
-{p[ ■ h +P?)V2{p't - kt)F{Q^)]f\\h) 

+ [{Vi{p't - kt) - 200oWt - kt)F{Q^)) 

-(m« +w')(p; • kt+pf)V 2 F{Q^)\fi{kt)) 

UJ 1 

1 — 2ujg{MA — ujujd — a;( — b/uP — Icos^pt) 
■v-p't^ {[{Viip't - kt) + 2uJDV2{p't - kt)F{Q‘^)) 
-{ms - uj'^){Vi{p't - kt) 

+^V2{p't-kt)F{Q^))]ff{kt) 

Pt 

+ [{ms - uj's){Vi{pt - kt) 


-2ojDV2{p't-kt)F{Q^))?!^ 

Pt 

-{v\ ■ kt+p?)V2{p't - kt)F{Q^)]f^{kt)}y ( 20 ) 


gi = 


h = g^ = = [ Fi 


Ma 

Mas 


92 = t2 = gs = ts = 


Ma, 


-F 2 


T 

93 


— 

to — 


F 2 

Ma, 


■{Mas-Ma), 


F 2 

Ma, 


'{Ma, + Ma), 


T 

9i 


= ti = 


F 


2 „2 


Ma, 


9 ■ 


(18) 


The BS wave function of Ab was given in previ¬ 
ous works and has the form Xp’’{p) = 4’^'’{p)ua,{v, s), 
where (f>^'’{p) is the scalar BS wave functions [171142] . 


F 2 



d?pt 

(27r)3 


^'^’’{Pt) 


dAkt 

(27r)3 

1 


1 — 2uj's{Ma — ojojd — ^'s — VoF-M cos 9Pt) 
« • Pj X I [{Vi{p't - kt) + 2uJDV2{p't - kt)F{Q^)) 


-{ms - w')(Ui(p; - kt) + ^U 2 F(Q 2 ))]/ 7 (fc,) 

Pt 

-^[(ms - u}'s){Vi{p't - kt) - 2wD%{p't - kt)F{Q'^)) 
^ - {p{ ■ kt+p?)V2{p’t - kt)F{Q^)]f^{kt)y 

( 21 ) 
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where uJs = yjml-pi, w' = yjrnl-p^^, p[ (= 
p'l ■ v) and p'l (= p' ■ v) are the transverse and 
tudinal relative momenta along the momentum 
respectively, oj = v ■ v' (v and v' are the velc 
of the Ah and A, respectively) is the invariant ve 
transfer, 0 is the angle between pt and . All the 
factors are functions of the invariant velocity tra 

uj = —^ - - -, therefore, the minimum and 

2mA^mA ’ ’ 

mum values of a; are 1 and 2 ^^ mA ’ I'^spective 
our calculation, we take = 0.45 GeV, Ma = 
GeV, Ma^ = 5.62 GeV. Then one can find w r 
from 1 to 2.62. Substituting Fi and F 2 into Eq 
we will get the numerical results of pi, ti, gf a: 
as functions of w. The plots of 51 , g 2 ^ g^, gj, a; 
are shown in Figs 1 and 2. Other form factors c 
obtained from Eq. (18) straightforwardly. 



Fig. 1 Ab ^ A form factors as functions of cj. The 
(lower) line in the upper plane corresponds to gi (gj 
upper, middle, and lower lines in the lower plane corr( 
to 32 , 9i', and tj, respectively. The solid and dotted lin 
respond to mn = 0.7 GeV and mn = 0.8 GeV, respec 
when K = O.OSGeV®. 


In a similar way, we obtain the form factors for 
p replacing rus by m„ and Ma by Mp. uj for Ab —f 
p ranges from 1 to 3.08, and the the Clebsch-Gorden 
coefficient of the u{ud)op configuration is l/-\/2- The 
numerical results for gi, 32 , tf , tj, and g^ for Ah ^ p 
are plotted in Fig. 3 and 4. 

From Figs. 1-4, we can see that the magnitudes of 
form factors decrease as ui increases. This is because the 
overlap integrals of BS wave functions decrease with the 
increase of uj. The numerical results of these functions 
will be used to calculate the decay widths of Ah —f 
Al^l~ and Ah —f plv in next sections. 



Fig. 2 Ab ^ A form factors as functions of uj. The upper 
(lower) line in the upper plane corresponds to gi (gj’). The 
upper, middle, and lower lines in the lower plane correspond 
to 32 , g^, and respectively. The solid and dotted lines 
correspond to k = 0.02 GeV® and 0.08 GeV®, respectively, 
when mo =0.75 GeV. 



Fig. 3 Ab ^ p form factors as functions of uj. The upper 
(lower) line in the upper plane corresponds to gi {g^). The 
upper, middle, and lower lines in the lower plane correspond 
to 32 , g^, and tj, respectively. The solid and dotted lines cor¬ 
respond to mo = 0.7 GeV and mo = 0.8 GeV, respectively, 
when K = 0.05 GeV ®. 


III. Ah —>■ Al~^l Decay in the Standard Model 

We first apply our results for the Ah ^ A form fac¬ 
tors to calculate the differeutial fraction for the decay 
Ah —f [l = e^pLjT) in the SM. This process is 

loop-suppressed, and hence potentially sensitive to new 
physics beyond the SM. At the quark level. Ah —f Al~^l~ 
is described by 5 —f transition. The effective 

Hamiltonian describing the electroweak penguin and 
weak box diagrams related to this transition is given 

by [nuns] 

+ Cw{s-i^PLb){ir^hl) 


GpOi 

TT 
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M{Ab -)> Al+l ) = X Xt[l'-ff^l{uA['y^{Ai{l + 75) 

2v 27 r 

+Bi{l - 75)) + ia^''pi,{A 2 {l + J5) 

+^2(1 - 75 ))]uyiJ 
+ hfj.j 5 l{uA[l^{Di{l + 75) 

+Ei{l - 75)) + ia^''p,j{D2{l +75) + £2(1 - 75)) 
+p^iD3{l + 75) + Esil - 75 ))]maJ] , ( 23 ) 

where the parameters Ai, Bi and Dj, Ej (i = 1,2 
and j = 1 , 2 , 3 ) are defined as 


Fig. 4 At, ^ p form factors as functions of ui. The upper 
(lower) line in the upper plane corresponds to gi (gj). The 
upper, middle, and lower lines in the lower plane correspond 
to 52 , di , and tj, respectively. The solid and dotted lines 
correspond to k, = 0.02 GeV® and 0.08 GeV®, respectively, 
when m_D = 0.75 GeV. 


—2C^^^mb 



ihpi) , 


( 22 ) 


where Gp ts the Fermi coupling constant, a is the 
electromagnetic coupling constant. 


Bi = {gi + ti) - (ffz^ -if)|, 

Dj = —Cio{gj — tj), 

Ej = -^ 10 ( 5 ^ +tj)- (24) 

The final task is to calculate the decay rate of Tl —>■ 
Al^l~ in the whole physical region, 4m^ < < (^dt ~ 

ttia)^- The differential decay rate is obtained as 


At = VtbV*, 

Pr,l = 2(1 ± 75 ), 

q is the momentum transferred to the lepton pair which 
is the sum of the momenta of and l~, , 

and Cio are the Wilson coefficients. We note that only 
the term associated with the Wilson coefficient Cio is 
independent of the renormalization scale. To find the 
transition amplitude A^ -A A, we need to sandwich 
this effective Hamiltonian between the initial and final 
baryon states and calculate the matrix elements 

(A(P',s')lsj^blAb(F,s)), 


(A(P', s')|s7g75&|4lb(f’, s)), 
(A(P',s')lsia'^-q^blAi,(P,s)}, 


and 

(A(P', s')lsia'^^j5q‘'blAi,(P, s)}. 

These matrix elements are expressed in terms of the 
form factors obtained from Eqs. (16-21) in Section II. 
Then, the matrix element of the decay Ai, —>■ Al~^l~ can 
be written as 


dP 

dq^ 


GW 

81927r^myi,^ 


IVtWfviWiGWlTois) 

+lr2is)], 


where 


A(l, r,s) 


s = q^/m\^, 
r = iriAlmAi, 

1 + r'^ + s'^ — 2r — 2s — 2rs, 


(25) 


and 



is the lepton velocity. The functions 7o(s) and 72 (s) are 
given as [laisg 


To{s) = 32mfmW{l + r- s)(|£> 3 p + l^^aP) 

-b64m?m^^(l - r - s)Re{DlEg + DW) 
+64m^^V^(6m;^ - M^^s)i?e(i7*£’i) 
+&‘lm'fm\\/r(2'mA,,sRe{D'^Eg) 

+(I - r + s)Re{DlD3 + EW)) 

-|-32m^(2mf -|- m\s){{l — r + s)mAtV^Re{AlA 2 
+BIB 2 ) - m^,(I - r - s)Re{AlB 2 + A*Bi) 

-2^{Re{AlBi) -(- m'^AsReiAW)) 
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+8m^j4m^(l + r- s) + + r)2 - 

s2)](|^iP + m^) + Sm^ j4mnA + (1 + r - s)s] 
+ml^s[{l-rf-s^]}(\A^\^ + \B^\^) 
-8m^j4m?(l + r - s) - myij(l - rf - s^]} 
(|£iiP + \Ei\'^) + 8m^^s7;2{-8mA,sVfiie(£i;£;2) 
+4(1 - r + s)^Re{DlD 2 + £;i£; 2 ) - 4(1 - r - s) 
i?e(i^r £;2 + D^Ei) + [(1 - rf - s^] 

(|i^2p + |i?2p)}, (26) 


, which is shown in Fig 5 along with recent ex¬ 
perimental results from LHCb [3] ■ The agreement of our 
results with the experimental data is clear except when 
the square of the momentum transfer is bigger than 15 
GeV^. So there is still room for possible new physics. 
We also compared our result with that of the first Lat¬ 
tice QCD simulation [HES] and found that they are 
consistent with each other. Predictions for Af, —>■ Al^l~ 
when I = e,T are shown in Fig. 6 and Fig. 7, respec¬ 
tively, which will be compared with the forthcoming 
experimental data. 


T2{s) = -8mluf A(|7lip + + |Cip + \D,\^) 

+8m^,s^?A(|7l2p + |i?2p + IC2P + 17521"). 

(27) 

oj = + m\ — q^)/2mA^rnA, so w ranges from 1 

to {niA + — 4mf)/2mAmA^- The differential decay 

rate expressed in terms of w has the following form, 


dr 

d{2mAmAt^) 


+ lr2{s)]. (28) 


In our numerical calculations, we use the value of the 
CKM matrix elements |VtbV)*| = 0.041 and the Wilson 
coefficients at /r = rn-b, = —0.313, = 4.334 

and Cio = —4.669 [SHESUSD]. As mentioned before, 
letting K range from 0.02 GeV to 0.08 GeV and m_D 
from 0.7 MeV to 0.8 MeV, we have obtained numerical 
results of the form factors g{t)i, {i = 1,2,3). 

Using the lifetime of the Ab baryon, (1.451 ± 0.013) x 
10“ [HI], and integrating the differential branching 
ratio (28) over a; from 1 to (m^ + m^^ —4TO^)/2myiTOyi,,, 
we obtain the ranges of the branching ratios, which are 
listed in Table I. 

In Table I, we also present the values of the branch¬ 
ing ratios obtained in HQET [54] and the light-cone 
QCD sum rules [18]. It can be seen from the Table I 
that, as is excepted, the branching ratios decrease when 
I goes from the e to r [HD] . We can also see that our re¬ 
sult on the branching ratio for Ab —>■ AgAgL~ is about 
I/4 of that predicted by HQET and about 1/6 of that 
given by light-cone QCD sum rules. Our result is consis¬ 
tent with the experimental data and the other two are 
not. Since 10^°^10^^ AbAb pairs are expected to be pro¬ 
duced per year at LHCb, the results presented in Table 
I indicate that the detection possibility of Ab —>■ Al^l~ 
{I = e,p,,r) is quite high [TS] . 

Letting k and mD vary in their regions we obtain 
the area of the differential branching ratio for Ab —>■ 



fiGeV^) 


Fig. 5 The differential branching ratio for Ab —>■ Ag~^ g 
obtained in our model. The upper (lower) curved solid line 
corresponds to the upper (lower) boundary of the differential 
branching ratio as mo and k vary in their ranges. The ex¬ 
perimental data are taken from Ref. [11], with the error bars 
including systematic and statistical uncertainties. 



Fig. 6 The differential branching ratio for Ab —>■ de+e 
obtained in our model. The upper (lower) curved solid line 
corresponds to the upper (lower) boundary of the differential 
branching ratio. The uncertainties of the boundaries come 
from the ranges of mo and re. 
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Table 1 Values of the branching ratios for A}, —>■ yll+1 in our model and the values from the light-cone QCD sum rules and 
HQET for different leptons 



present work 


HQET [54] 

light-cone QCD sum rulesll8| 

Exp. [51] 

Br{Ab -A yle+e ) (1.21 

- 2.32) X 10“® 

(2.23 - 

- 3.34) X 10-® 

(4.6 ± 1.6) X 10-® 


Br{Ab ^ Ap,+p,~) (0.53 

- 0.89) X 10-® 

(2.08 - 

- 3.19) X 10-® 

(4.0 ± 1.2) X 10-® 

(1.08 ±0.28) X 10- 

BriAb^Ar+T-) (0.037 - 

- 0.083) X 10-® 

(0.179 - 

0.276) X 10“® 

(0.8 ±0.3) X 10-® 




Fig. 7 The differential branching ratio for ylt —^ ylr+r 
obtained in our model. The upper (lower) curved solid line 
corresponds to the upper (lower) boundary of the differential 
branching ratio. The uncertainties of the boundaries come 
from the ranges of mo and k. 


IV. THE DECAY At -A- pW 

At the quark level, the decay Ab ^ plv is governed by 
the tree-level b ^ u transition. The effective Hamilto¬ 
nian responsible for this transition at the quark level 
has the form 

Ti{Ab pW) = ^VubUAfii^ - A5)bh^{_l - (29) 

To calculate the amplitude, we need to sandwich the 
above Hamiltonian between the initial and final states 
and compute the matrix element (p|m 7^(1 — ^b)b\Ab)- 
As discussed in Sect. I, we have obtained the six form 
factors for this decay amplitude in our previous work 

m- 

In the next step, we calculate the total decay rate 
of Ab —f plv in the whole physical region, 

mf <q^ < - irip)^. 

The range of w for this decay is 

[1, {m% +m\- mf)/{2mA,mA)] 


The decay width is given by the following expression; 
r{Ab^pW) = 384^3^^3 “ rnf/qY 

X^{S^-q^){A^-q^)N{q% (30) 

where 

N{ q^) = glq\A\Aq^ - mf) + A\l + 2m^/q^) 

+ 2q^){2q^ + mf )) + gliq^KA^ - q^) 

(2S^ + q^) X {2q^ + mf) + 3Ml^gl{q^)m^{S^ 

-q^)q^ + 6g,iq^)g2iq^)iA^ - q^){2q^ + mf)S 
- 6 ffi(g")ff 39 "mf (A^ - q^)A + t\{q^){S\Aq^ - mf) 
+2E^A'^{1 + 2ml)/q^ - (A^ + 2q^)(2q^ + mf)) 
+tl(q^)(i:^-q^)(2A^ + q^)(2q^ + mf) 
+3tl(q^)mf(A^ - q^)q^ - %h{q^)t^{q^){E^ - q^) 
{2q^ + mf)A + 6 ti(g 2)53 + {q'^)mf{A'^ - q^)E, (31) 

and A = niA^ — nrip, E = mA^, + rup. The numerical 
results are listed in Table H, together with theoretical 
results from other models. 

The four results of Ref. [56] in each line in the table 
refer to those from QCD sum rules, lattice QCD, QCD 
sum rules in the heavy quark limit, and lattice QCD in 
the heavy quark limit in order. The decay rates were 
calculated in the covariant quark model [ST], SU{3) 
symmetry quark model [35] and HONR and HOSR con¬ 
stituent quark models (HONR and HOSR refer to har¬ 
monic oscillator nonrelativistic and harmonic oscillator 
semirelativistic constituent quark models, respectively) 
[59] . We compare our results for the rates [in the units of 
\Vub\^ ps~^] with the predictions of other phenomeno¬ 
logical methods. From the table, it is clear that our 
results are of the same order as those of lattice QCD 
in the heavy quark limit [56] and those from Refs. [58- 
59]. However, our results disagree (up to two orders of 
magnitude) with those of Ref. [57] and those obtained 
from QCD sum rules, lattice QCD, and QCD sum rules 
in the heavy quark limit [SSj. 

V. SUMMARY AND DISCUSSION 

Theoretical studies of the rare baryon decay of Ab -A 
Al^l~ require knowledge of the hadronic matrix ele- 
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Table 2 Values of the decay rates (in units | ps of the transitions At —>■ plu and comparision with other methods 


present work Other models 

At ppi-vt, 3.12 ~ 7.06 250 ± 85[5S];235 ± 85[5S];477 ± 175(55]; 3.84 ± 1.25(55]; 13.3(57]; 6.481551:4.551551:7.551551: 

At pe-v^ 2.09 ~ 7.40 6.48(5S];4.55(5n];7.55(51]; 13.3(57]; 250 ± 85(55]; 235 ± 85(55]; 478 ± 175(^; 3.76 ± 1.20(55]; 

At pr-p^ 1.35 ~ 4.10 4.01(51]; 6.55(51]; 312 ± 105(55]; 208 ± 70(55]; 646 ± 215(55]; 1.93 ± 0.70(55]; 9.61571: 


ment {A\srh\At) which involves nonperturbative QCD 
effects. At the leading order in HQET, this matrix el¬ 
ement is described by two independent form factors, 
which are determinated by the wave functions of the ini¬ 
tial and final baryons. We calculate these two form fac¬ 
tors in the BS equation approach in the quark-diquark 
model. Consequently, we obtain all the twelve form fac¬ 
tors resposible for the decay At -A Al'^l~ , which de¬ 
pend on the two parameters, niD and k, in our model. 
Then, we obtain the total and the differential branch¬ 
ing ratios of the decay At -A Al'^l~. We also compare 
our results with those of other approaches and the ex¬ 
perimental data from LHC. We find that our result on 
the total branching ratio of At -A ApL^pL~ is consistent 
with the experimental data but those of the light-cone 
QCD sum rules and HQET methods are not. We also 
obtain the area of the differential branching ratio for 
At -A ApAiJL~, which is consistent with those of the 
first lattice QCD simulation and with the experimental 
data except when the square of momentum transfer is 
bigger than 15 GeV^. This indicates there is still room 
for possible new physics effects. Eurthermore, we give 
the decay rates of At -A plD. We find that our values 
are the same order as those of Refs. [58-59] and those 
of lattice QCD in the heavy quark limit in each line 
in Table H (^ , but different from those obtained from 
QCD sum rules, lattice QCD, and QCD sum rules in 
the heavy quark limit by up to two orders of magnitudes 
[5S] . The decay At -A plv will likely yield the first de¬ 
termination of the CKM matrix elment \Vut\ from At 
decays at LHC. Our results depend on two parameters 
in our model, mu and k, which vary in some ranges. 
This lead to some uncertainties in our results. All our 
predictions will be tested in the future experiments. 
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